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Abstract 


An atmospheric flight vehicle in hover is typically controlled by varying its thrust 
vector. Achieving both levitation and attitude control with the propulsion system 
places considerable demands on it for agility and precision, particularly if the vehicle 
is statically unstable, or nearly so. These demands can be relaxed by introducing 
an appropriately sized angular momentum bias aligned with the vehicle’s yaw axis, 
thus providing an additional margin of attitude stability about the roll and pitch 
axes. This paper describes a methodical approach for trading off angular momen- 
tum bias level needed with desired levels of vehicle response due to the design 
disturbance environment given a vehicle’s physical parameters. It also describes 
several simplifications that provide a more physical and intuitive understanding of 
dual-spin dynamics for hovering atmospheric vehicles. This approach also mitigates 
the need for control torques and inadvertent actuator saturation difficulties in try- 
ing to stabilize a vehicle via control torques produced by unsteady aerodynamics, 
thrust vectoring, and unsteady throttling. Simulation results, based on a subscale 
laboratory test flying platform, demonstrate significant improvements in the atti- 
tude control robustness of the vehicle with respect to both wind disturbances and 
off-center of gravity payload changes during flight. 


1 



1 Introduction 


There has been recent interest in thrust-levitating flight vehicles capable of vertical takeoff 
and landing (VTOL) for reconaissance [1], [2], [3], [4], and more exotic applications such as 
rescue platforms for tall buildings [5], air ambulances [6], flying cars [7], and an air scooter [8]. 
These vehicles typically control their attitude in hover by varying their thrust vector. This 
requires agility and precision in controlling the thrust vector and, additionally, an adequate 
mathematical model of the thrust vector dynamics for designing the control law. 

In spacecraft, propulsive attitude control consumes scarce expendables and can contaminate 
the spacecraft’s instruments. Therefore, dual-spin stabilization is frequently preferred. Dual- 
spin refers to a partitioning of a body into a component that spins about an axis that is 
fixed in the body, and a component that remains stationary in body axes. The presence of 
the spinning component results in an angular momentum vector of fixed orientation for the 
entire system. This angular momentum bias provides a restoring moment against attitude 
perturbations orthogonal to the spin axis. Such passive stabilization is identical to that seen in 
spinning toy tops. For more detail description and discussion on dual-spin systems for spacecraft 
applications, namely, “gyrostat” or “bias momentum satellites”, see for example, Section 6.6 in 
[9], Chapter 5 in [10] or Chapter 11 in [11], 

For thrust-levitated vehicles, the objection to relying entirely on thrust vectoring for attitude 
stabilization primarily centers on the difficulty of modelling the response of the thrust vector 
to commands, because of unsteady fluidic phenomena. This difficulty is exacerbated when 
the vehicle is statically neutrally stable or unstable, since the thrust variations required for 
stabilization are now more active. Therefore an accurate model is needed for control law design 
in the presence of unsteady flow phenomena during thrust vectoring. Inclusion of a dual- 
spin mechanism in such vehicles, operating along the z-body axis (down), can ameliorate this 
difficulty by passively enhancing the stability of the vehicle, reducing the demands on the thrust 
vectoring system. 

The disturbance environment and maneuvering envelope of thrust-levitating air vehicles 
is significantly different from that of spacecraft and therefore the criteria for specifying the 
minimum size of the dual-spin system’s angular momentum bias are also different. For dual- 
spin spacecraft, a design practice as discussed in Section 15.3 of [9] for sizing the level of angular 
momentum is based on the anticipated change in the angular momentum in the system due to 
unknown external torques, i.e., 



ft 

1 » 

/ T d ,st dt 

Jto 


If the disturbance torques are cyclic in nature and the resulting perturbations are less than 
mission requirements, then the above passive stabilization will be adequate. However, if the 
disturbance torques have a non zero mean or biased trend (i.e. secular disturbances) and 
exceed mission attitude constraints, then active control will be additionally needed. This paper 
examines in detail the basic requirements and sizing of bias momentum for a generic thrust- 
levitating vehicle system, and the related aspects of the dynamics that govern a dual-spin 
system. 
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Section 2 briefly discuss relevant dynamics for a generic dual-spin system and introduce a 
gyroscopically dominant state that leads to analytical simplifications without explicit lineariza- 
tion assumptions. Section 3 examine the pitch and roll response characteristics due to pitch 
and roll disturances for a dual-spin system about a hovering trim condition. In particular, the 
precession frequency and damping for a generic dual-spin system are defined and the effect 
of the bias momentum level on the system response to disturbances is examined. Section 4 
gives a derivation of an approximate analytical expression for the worst case response bound 
using maximum singular value of the frequency response matrix. Its dependence on the levels 
of disturbance, precession damping, and bias momentum are clearly shown. In Section 5, we 
examine mean square response to disturbances quantified by power spectral densities. As in 
Section 4, approximate analytical expressions are derived for the mean square response under 
various conditions. Section 6 focuses on the bias momentum sizing issue for a generic flying test 
platform and examines the dependence of the vehicle response on various parameters, and pro- 
vides rationale for a recommended design range for bias momentum levels. Section 7 presents 
simulation results of the generic dual-spin vehicle under simulated wind disturbances and in- 
ertia changes during flight. The simulation includes a validation of the approximate response 
formulae used for bias momentum sizing. The simulation results show significant improvements 
in the robustness in the roll and pitch responses to large wind disturbances and inertia changes. 
Conclusions are given in Section 8. 
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2 Dynamics of a dual-spin system 


Consider a dual-spin system TZ + B which consists of a main rigid body platform 7 Z, and a 
flywheel B , as illustrated in Figure 1. For simplicity, let TZ , B , and TZ + B, be centered about 


A 



Figure 1: Schematic of a dual-spin system. 

a common nominal center of mass, O cm . Furthermore, using the vectrix notation in [11], let 
F t — (bi, b 3 , b 3 ) denote the unit triad corresponding to the principal inertia axes for TZ , B, and 
TZ + B. In the following paragraphs, the use of the vectrix notation allows one to succiently 
assemble matrices of vectors and perform mathematical operations including vector operations 
and matrix algebra. 

Let cu denote the inertial angular velocity of TZ. Suppose the flywheel, B , is axi-symmetric 
with inertia dyadic, / and spins at a rate VLb with respect to the main platform, TZ , with 

inertia / , along its b 3 axis which is fixed in TZ. Hence, the inertial angular velocity of B, can 

be written as uj = u + QBb 3 . The angular momentum of 7 Z + B about its center of mass can 
then be written as 

h — I ■ uj + I ■ uj (2.1) 

— > — >R —> —>B — >B v ’ 

By defining the net external torque disturbance on TZ + B about its center of mass as 


_r = F 1 r, t = { T 2 
t 3 

the rotational equations of motion for TZ + B in vector form 


— — ( h ) + u x h 


( 2 . 2 ) 


( 2 . 3 ) 
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can be expressed in the following scalar form 


T\ 

= hio i + uj 2 [h B + a ; 3 (/ 3 - I 2 )\ 

(2.4) 


= I 2 u 2 — uji [h B + u ; 3 (/ 3 — ii)] 

(2.5) 

r-s 

= I 3 UJ 3 + h B + ( I 2 — I\)u\U 2 

(2.6) 


where the following inertia dyadics and vectors are expressed in principal inertia coordinates, 
as follows: 


I 4 1+1 = F t IF, where I 4 diag (/ ls I 2 , J 3 ) 

* ' Jtx ^ Jd 

4 b — F T IbF, where I B = diag (i^, 1^, I a B ) 

4 R - fTi R F \ where/* = diag (l l R , 1% 1%) 

A r -i l 

cu = F t o;, where u; = < u ;2 > 



where h B = I R fl B 


hi + Ib 


(2.7) 

(2.8) 
(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


Note that if the angular momentum of flywheel due to its spin, h B , is zero, then the above 
dual-spin dynamical equations 2.4 to 2.6 reduces to the familiar Euler equations describing the 
attitude dynamics for a single rigid body (see for example [12], [13]). The terms, u 2 h B and 
uqhg, appearing in the roll (1-axis) and pitch (2-axis) equations 2.4 and 2.5, represent the 
internal gyroscopic torques due to the bias momentum in the flywheel. Clearly, these internal 
gyroscopic torques can be made to “dominate” the motion with a sufficiently large angular 
momentum in the dual-spin system. Notice in equation 2.6 that a similar bias momentum 
induced gyroscopic torque does not appear in the yaw axis, rather, the reaction torque, h Bl 
due to an angular acceleration of the flywheel appears. It is well known that the presence (and 
absence) of the bias momentum induced gyroscopic torques in the roll and pitch axes (in yaw 
axis) gives rise to the directional stability (unrestrained motion) in the roll and pitch axes. 


Suppose the flywheel spins at a constant rate 

\h B \ = 0 (2.12) 

and the bias momentum in the flywheel dominates in the following sense: 


Definition 1 The dual-spin dynamical system defined by equations 2-4 to 2.6 is said to be a 
“gyroscopically dominant ” state if 

\h B | » max (|(/ 3 - J 2 )|, |(/ 3 - h)\) ■ M (2.13) 


Then, the dynamics of this generic dual spin system can be approximated by 

fidJi « -c o 2 h B + T\ 

I 2 io 2 Rj U\h B + r 2 

= ~ (I 2 ~ Il)cOiU 2 + 73 
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(2.14) 

(2.15) 

(2.16) 


Alternately, if the dual spin equations in 2.4 to 2.6 are linearized about angular velocities of 
Ui = 0 by assuming \ui\ -C 1, for i = 1, 2, 3, the following results 


IlUJl 

^ -U 2 h B + Ti 

(2.17) 

1-2^2 

~ Wi h B + t 2 

(2.18) 


~ t 3 

(2.19) 


Notice that the roll and pitch dynamics arrived at by bias momentum dominance, as given by 
equations 2.14 and 2.15, are identical to the corresponding linearized equations 2.17 and 2.18. 
However, the yaw dynamics are different in that equation 2.16 is an exact nonlinear equation 
while equation 2.19 is an approximate linear equation. 

In summary, as long as a dual-spin system is gyroscopically dominant , i.e., equation 2.13 
is valid, the linear roll and pitch dynamics, as given by equations 2.14 to 2.16, will be valid 
without having to assume small “linear” angular velocity motion. On the other hand, if the 
bias momentum level in the dual-spin system is such that it is not gyroscopically dominant, 
then linearity assumptions must be made in order for linear equations 2.14 and 2.15 to hold. 


6 


3 Pitch-Roll Dynamics about Trim 


As evident in equations 2.14 and 2.15, the roll and pitch states, oq and a ; 2 , are independent of 
the yaw state, uq, although the pitch and roll nonlinear gyroscopic term, (J 2 — ii)uqu; 2 influences 
yaw motion directly. With this fact, consider only the following linear roll and pitch dynamics 
from equations 2.14 and 2.15 


where 


A 


t J I ^2 J 


Cl 

h B 

r 

,h 

h B 

h 

C2 


. h 

h . 

1- 



A' 1 

0 



(3.1) 


(3.2) 


Notice from the diagonal elements of the system matrix A in equation 3.2 that the viscous 
damping terms <q and C 2 have been added. These terms are intended to represent the physical 
energy dissipation mechanism during the precession/nutation motion which are described by 
their respective angular velocities of uq, and u; 2 . For example, the energy dissipation associated 
with a “wobble” or “nutation” damper (see for example [14]), whether implemented passively 
or actively, can be accounted using the above viscous damper model. Although in general, 
damping or energy dissipation effects for a particular physical system are difficult to model 
accurately, in this study, these parameters are used in the limited context of a performance 
trade study for the purpose of momentum wheel sizing. 


In transfer function form, the above roll and pitch dynamics can be expressed as 


- G(s ; h ) 

W2 l r 2 


where 


G{s-h B ) = 


A (a) 


_S_ 1 C2 __ _B_ 

h "T hh I1I2 

JL 1 Cl 

hh h ^ hh J 


; A( S ) 4( a -A)(a-A*) 


with the complex conjugate eigenvalue pair 

A, A* = - < 2 


(3.3) 


(3.4) 


(3.5) 


and 



The transfer function matrix G(s]h B ) defines the roll and pitch angular velocity frequency 
response to external roll-pitch torque disturbances. The complex conjugate pair of eignevalues 
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in equation 3.5 corresponds to the often dominant precession motion due to the bias momentum 
induced gyroscopic coupling in the roll and pitch axes. This eigenvalue pair is parameterized 
in the familiar form for a second order oscillatory system involving damping ratio, (, and 
frequency, u 0 , given by equations 3.6 and 3.7, respectively. The term A 0 as given in equation 
3.8 denotes the undamped precession or wobbling frequency of oscillation and is analogous to 
the precession frequency in a dual-spin spacecraft as described for example in Chapter 6 of [11]. 
Notice from equation 3.7 that the damped frequency, u a , reduces to this precession frequency, 
A 0 , in the absence of damping, i.e., ( = 0. This precession frequency increases linearly with the 
level of bias momentum and is inversely proportional to the rotational inertia of the vehicle. 

Defining the following special forms for the damping ratio, (, as given in equation 3.6, will 
lend clarity in the following paragraphs: 



— if Ci — c 2 — c, I 1 = / 2 = /, and \c/h B \ <C 1 (3.11) 

B 


Notice from equation 3.11 that for “small” damping, the ratio, C m = can be assumed to be 
the damping ratio in the system. 

Figure 2 shows the log magnitudes of the roll and pitch angular velocity responses to corre- 
sponding disturbance torques, at bias momentum levels of (0, .1, 1, 10, 50) Nms. Small damping 
constants corresponding to C\ = C 2 = 0.02 h B are assumed in the above figures, corresponding 
to about 2 percent damping ratios for each case. The case with h B = 0 is shown by dashed 
lines in the figure. Note that there are no dashed lines in the off-diagonal plots since there is 
no coupling in roll and pitch. In all the subplots in Figure 2, notice that the responses decrease 
and their peaks shift to higher frequencies, with increasing bias momentum levels. Notice that 
increasing the level of bias momentum can significantly reduce the roll and pitch angular ve- 
locity response to disturbances. Furthermore, as the level of bias momentum increases, the 
roll response to roll disturbances decreases more rapidly than the roll response to pitch distur- 
bances. Alternately, the dynamically uncoupled system at zero bias momentum transforms to 
a system whose characteristic dynamics are gyroscopically dominant (see definition in equation 
2.13). 

Based on equation 3.3, Figure 3 shows the maximum singular value frequency response of 
the roll and pitch angular velocities, and is given as follows 

a [G(ju; h B )\ = sup Vi/ (3.12) 


where r = (ti,t 2 ), u = (uq,u; 2 ), and || • || 2 denotes the Euclidean norm on T n . If the above 






Figure 2: Roll and pitch angular velocity responses to roll and pitch disturbance torques for 
bias momentum levels of (0, .1, 1, 10, 50) Nms; c\/h B = C 2 /h B = 0.02. 
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supremum of Eulidean norm ratio is searched for the worst case over all frequencies, 


maxcr [G(jw, h B )\ = || G(jv; h B ) ||oo (3.13) 

V 

then the above norm can also be expressed in terms of induced signal 2-norms (see for 
example [15], [16]). The maximum singular value captures the directionally worst case response 



Figure 3: Maximum singular value angular velocity responses to roll and pitch disturbance 

torques for bias momentum levels of (0, .1, 1, 10, 50) Nms; C\/h B = c 2 /h fl = 0.02. 

due to simultaneous disturbance torques in roll and pitch. In both figures, the dashed line 
represents the reference case with near zero bias momentum. As in the previous plots, the 
responses shown in Figure 3 decrease and their peaks shift to higher frequencies, with increasing 
bias momentum levels. Notice from the maximum singular value frequency response that it 
looks similar to an undamped second order system with a resonance peak at the precession 
frequency given by equation 3.5. 

The next two sections consider two bias momentum sizing approaches for the generic dual- 
spin system which are based on the angular velocity responses of the main platform to a set of 
disturbances. These approaches differ in their assumptions in the disturbances and the metric 
used to describe the roll-pitch responses. As will be evident, various simplifying assumptions 
are made to obtain corresponding intuitive analytical expressions for responses which are only 
approximate but nevertheless useful for sizing and conceptual design of effective bias momentum 
levels. 
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4 Worst Case Response Bound 

4.1 Disturbance Model 


In this approach to sizing, assume that the external disturbances are unknown but bounded 
and can be expressed as 

T = V(ju)w (4.1) 

where V(ju) € C 2x2 denotes a frequency response matrix representation of a filter, which is 
driven by an unknown fictitious source signal vector, w € C 2xl , which is assumed to be unity 
bounded as follows 

\\w(ju)\\ 2 < 1, Mv (4.2) 

Suppose Ti and r 2 are independent disturbance signals having a flat frequency response over a 
bandwidth, (BW) T . Hence, from equations 4.1 and 4.2, the magnitude of the disturbance filter 
can be approximated as follows 


I Vik (jv) I 


h(jv)| 

Wkijv) I 


( 2 ^ 

0 , 


i — k 
i 7 ^ k 


Vz/ < (BW) T , i, k 


1,2 


(4.3) 


and the disturbance filter is assumed to roll-off rapidly outside the bandwidth. 


4.2 Static Approximation 


The frequency response matrix from equations 3.4 and 4.3 can be written as 


G{jv\h B )V{jv) 


1 

A(jV) 


(g + ^n 


hh' 

h B T/ 

hh' 11 


h B T/ 

'TJ2V22 




22 


(4.4) 


where A(ju) denotes the value of A(s) as defined in equation 3.4 but evaluated at s = jV. 
Based on the maximum singular value frequency response given by equation 3.12 and plotted 
in Figure 3, the following approximations are considered 


a[G(ju-h B )] « a[G(0-h B )\ , for z/ < A 0 (4.5) 

Since V(ju) is assumed to be diagonal with constant frequency responses as described in equa- 
tion 4.3 over a certain bandwidth, (BW) T , it follows that over frequency, u, where 

u:u< {BW) r < A 0 (4.6) 
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the maximum singular value frequency response can be approximated by the following static 
approximations: 


a[G(ju-h B )V(ju)] 


a[G(0-h B )V(0)] 
1 


M 1 + W 


a 


f-Vn - V 2 


22 




max 


Wi + ff 


L> + v 7 H 2 - E, yD — y/D 2 — E) (4.8) 


max ( | Vii | , | V 22 I), if ci = c 2 = c 


k B\l 1 + ( 


(4.7) 


(4.9) 


where 


D = 


1 

2 


!+(£ 


1 1 C\C2 

1 I h 2 




Vii + x 


! + (£ 


! + lf 




22 


E = VfiVzl 


(4.10) 

(4.11) 


Notice that with the additional assumption in the damping constants satisfying the relations 
C 2 = Ci = c, the response formula simplifies considerably to equation 4.9. Furthermore, if the 
system is undamped or is only lightly damped, i.e., 


C\/h B -C 1 and c 2 /h s <C 1 or c/h B -C 1 (4-12) 

then the damping term can be completely ignored, simplifying even further. 

Consider bounding the maximum frequency response of the roll and pitch angular velocities 
given a bounded set of disturbance torques as described in equations 4.2 and 4.3. The roll-pitch 
angular velocity response from equations 3.3 and 4.1 is 


u(jv) = G{ju ; h B )V(ju)w(ju) (4.13) 

Based on the worst case response at each frequency, as given in equation 4.9, quantify the above 
frequency response bound requirement approximately as follows 

sup |KjV)|| 2 = a [G(ju; h B )V{jv)] « (l^nU < || oj target (u)\\ 2 (4.14) 

The term ||oy arsei (z/) || 2 denotes a, maximum desirable angular velocity response at the frequency 
of interest, u. 
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5 Mean Square Response Bound 


In this alternate approach, consider the mean square response (MSR) in the roll and pitch 
angular velocities given that the disturbance torques are random processes with known power 
spectral densities (PSD). 


5.1 Disturbance Model 


The disturbance filter model in equation 4.1 is assumed except that a different meaning is 
attributed to the fictitious signal w, assumed now to be a unit white noise process so that its 
PSD for i — 1,2, is 

Swiiy) = 1, Vi/ (5.1) 

As in the previous approach, assume that the disturbance torques, Ti, and r 2 are independent 
and have a near constant spectral density over a certain bandwidth, (BW) r . This means that 
the disturbance filter is diagonal and has the same form as described in equation 4.3. 

With the above unit white noise driven disturbance torque model, since their PSD’s are 
related by 

S Tl {v) = mjv)\ 2 S^v) (5.2) 

the frequency response magnitude corresponding to equation 4.3 is 


I V ik (jv)\ 


hO)l 

Wkijv) I 


y/SM, 

0, 


i — k 
i 7 ^ k 


Vi/ < (BW) t , i,k 


1,2 


(5.3) 


Analogously, assume that the disturbance torque PSD rolls-off rapidly beyond a certain band- 
width, (BW) t . 

Suppose the PSD of the anticipated disturbance torques, S Tl ( jv ) and S T2 (ju), can be approx- 
imated by constants, S° and S° 2 , over a bandwidth, BW r , as given in equation 4.6 and rolls-off 
rapidly beyond, BW r . For instance, in dealing with gust loads on aircraft, Von Karman and 
Dryden gust velocity PSD’s are used widely in design and simulation. These spectra basically 
approximate flat, low pass filters whose parameters depend on various flight conditions such as 
nominal airspeed and altitude. The break frequencies for these well known gust PSD’s are low 
(for example .1 to 1 Hertz) and they roll-off at orders 5/3 and 2 respectively (see for example 
[17]). These assumptions on the gust PSD structure lead to the following simplifying formulae 
for the mean square torque 

£[t t t]»2(S« + S° t1 )(BW) t (5.4) 
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5.2 Static Approximation 


For the system given in equation 4.13, the response PSD for the roll-pitch angular velocities 
can be expressed as [18] 


S u (”) = 


SM 

SM 


H n | 2 \h 12 \ 2 i r S^u) 
H 21 \ 2 \H 22 \ 2 \ \ S^{y) 


(5.5) 


where 

Hik(jv) = [G{jv\ h B )V(ju)\ ik (5.6) 

denotes the (i, k) frequency response from fcth-white noise disturbance torque to ith-angular 
velocity of vehicle. The mean square response of the zth-angular velocity response can be 
obtained from the response spectra given in equation 5.5 as 



(5.7) 


and with the white noise in equation 5.1, the disturbance spectra given in equation 5.3, and 
the coupled roll-pitch vehicle dynamics given in equation 3.4, the above mean square response 
expression, after a little algebra, takes the form 





1 

IaMF 

i 

I a O)I 2 



The mean square response of the roll-pitch angular velocities can then be written as 


(5.8) 

(5.9) 



(5.10) 

(5.11) 

(5.12) 


For sizing purposes, consider simple yet physically consistent approximations to the mean 
square formula given in equation 5.10. In addition to assumptions in equations 4.6 and 5.4, 
by limiting the disturbance bandwidth with respect to the precession frequency as given in 
equation 4.6, i.e. , 


- < 1 (5.13) 
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terms with the frequency ratio, A-, appearing in the integrand in equation 5.10 can be ignored. 

Ao 1 

In summary, if one approximates the integrand in equation 5.10 by its value at v — 0, i.e. , its 
frequency response at low frequencies, one can readily obtain 


£ [uj t lu\ ~ 


h 2 

B 


! + ff 



( 

i + (A ) 2 

£{rl} + 

i + (A 2 

V 

L k B \ 


L K \ 


£[t t t ] if ci = C 2 = c 




if ci = C 2 = c and small damping ( < 1 


(5.14) 

(5.15) 

(5.16) 


Notice that the approximation formulas in equations 5.14 to 5.16 can be obtained directly from 
the state equations in 3.1 without considering the integral given in equation 5.10. This is done 
by (1) neglecting inertial torques/angular acceleration terms to obtain the approximate formula 
in equation 5.14, (2) additionally assuming ci = C 2 = c to obtain equation 5.15, and (3) addi- 
tionally assume small damping to obtain equation 5.16. All the above approximations assume 
a low frequency response wherein inertial torque or angular acceleration effects are sufficiently 
small such that the angular velocity response is dominated by the gyroscopic “stiffness” effects 
and precession damping. Analogous to a spring-mass-dashpot system, notice that the the gy- 
roscopic “stiffness” appear as conservative restoring torques, but they are proportional to cross 
angular velocities, lu±, and 1 x 2 , and the level of bias momentum is analogous to the degree of 
“stiffness”. The approximate formulae in equations 5.14 to 5.16 do not depend on the inertia 
of the vehicle, A, and / 2 , because the inertial effects are assumed small at low frequencies. 
Only the gyroscopic “stiffness” term, h B , approximately defines the angular velocity response 
to torque at low frequencies. 


5.3 Explicit Integration Approximation 

The following sections derive and examine response formulae from explicit integrations for two 
classes of systems: (1) with no damping, and 2) viscously damped. 


5.3.1 Undamped System 


For the zero viscous damping case, the MSR of the angular velocities 
written as 


1 


e[u u ;} = — 


K Joe pU v ) 


S n +S n + [^S n + ^S n ](^-) 2 )du, 
1 2 


where 


x 0 


BW r 

77 


from equation 5.10 is 
x 0 < 1 (5.17) 

(5.18) 
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Under the assumption that the PSD of the anticipated disturbance torque is approximated by 
constants, S ^ and S° 2 (see discussions surrounding equation 4.6), the integral of equation 5.17 
is solved using indefinite integrals as described in the Appendix and is written as 

e[uj t uj\ = ^-f(ki + k 2 ) 2 + 0-5(fci - k 2 ) ln( j ^ X ° ) V x Q < 1 (5.19) 

h B \ 1 x a x a J 

where 

h =S° n +S° n , h±TS° n + yS° n (5.20) 

U 1 2 

Consistent with physical expectations where an undamped mode excited harmonically will re- 
sult in an instability, the condition x D = 1 in equation 5.19 predicts an infinite response in MSR. 
Clearly, for an undamped system, the disturbance bandwidth cannot equal or exceed the pre- 
cession frequency so that the undamped response integral, equation 5.17, and its approximate 
solution, equation 5.19, are valid only for BW T < X a . 


Special Case A = / 2 = / 

to 


In this case, equation 5.19 with the aid of equation 5.4 reduces 





(5.21) 


5.3.2 System with Damping 


For the non-zero viscous damping case, the integral expression in equation 5.10 is solved ex- 
plicitly using indefinite integrals outlined in Appendix and is written as 


£\u u\ = 


Ao 

K 


k?, 


4 /? va +/? 2 

h 


- k , In 


2Cm V09 + /3 2 


+ h 


(xo + fl ) 2 + C; 
(x 0 - (3) 2 + Cr 
-1 (x 0 - (3 


tan 


Cl 


+ tan 


-i 


x 0 + P 

Cm 


where (3, ( m , and k? are defined in equations 5.12, 3.6, and 5.20 respectively and 


(5.22) 


^3 — [1 + (c 2 /h B ) 2 ]S° + [1 + (ci//z s ) 2 ]SU 


(5.23) 


Since the explicit formula in equation 5.22 is nice but difficult to visualize, consider the following 
simplifying assumptions: I± = J 2 = / and C\ = c 2 = c in order to gain more physical insight. 
In section 6.1.1, estimates of these errors are computed by numerical simulation and it appears 
that even if these assumptions are not exactly satisfied, these approximations appear to be 
sufficiently accurate for design and sizing purposes. With these assumptions the following 
parameters simplify to 


h = [1 + (c/h g ) 2 )k 2 , Cm = c/h B , f} = l 


(5.24) 


so that after substituting for the mean square torque in equation 5.4, equation 5.22 reduces as 
follows: 


£[lu t lu] 


2cI(BW) 1 


tan 


-i 


x 0 -l 


+ tan 


-l 


x 0 + 1 

Cm 


(5.25) 
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Special Case ( m C 1, i„ < 1 : This case applies for small damping and where the distur- 

bance bandwidth is less than its precession frequency. For this case, the response in equation 
5.25 simplifies to 





(5.26) 


The formula predicts that MSR varies with the inverse square of the bias momentum level. As 
the disturbance bandwidth approaches the lightly damped precession frequency from below, the 
MSR is predicted to increase by the inverse factor of {l — x 2 0 ). With small disturbance bandwidth 
and large precession frequency or bias momentum, the above formula for MSR matches both 
the static approximation derived previously in equation 5.16, and the undamped special case 
based on explicit integration approximation given in equation 5.21. 


Special Case ( m -C 1, x Q > 1 : This case applies also for small damping but where the 

disturbance bandwidth is greater than its precession frequency. For this case, equation 5.25 
simpifies to 


£[uj T (jj\ 


7T 


2cI(BW) 




(5.27) 


When disturbance bandwidth is larger than precession frequency, the formula predicts that 
the MSR of the angular velocities will not be dependent on bias angular momentum. This is 
because with small damping, the MSR will be dominated by the large response in the lightly 
damped precession mode. As long as the disturbance bandwidth equals or exceeds the preces- 
sion frequency, the contributions to MSR from higher frequencies will be small. The formula 
also predicts that larger inertia platform and damping levels will result in smaller MSR to 
disturbances. Finally, the MSR is also inversely proportional to disturbance bandwidth which 
reflects the fact that for a given mean square disturbance torque, its spectral density, inversely 
proportional to its bandwidth, will be spread out well beyond the bandwidth of the system, so 
that the system response will be less with increasing bandwidth of the disturbance. 


Special Case 1, any x Q : This case applies for large damping for any disturbance 

bandwidth. Using the Taylor expansion of tan -1 written in Appendix, equation 5.25 can be 
approximated as 


£[uj t uj] 


e[r 1 


t 


(5.28) 


With large damping, equation 5.28 predicts that the MSR of the angular velocity will be 
dependent only on the damping coefficient. Physically, this is the case where the damping 
torques in the system dominates over both gyroscopic and inertia torques. The analogy to a 
spring-mass-dashpot system is a response to disturbances when both the mass and stiffness in 
the system are very small, a very unlikely scenario in the scope of applications of interest in 
this study. 
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6 Bias Momentum Sizing 


In general, the sizing of bias momentum to improve directional stability depends on many factors 
and parameters including vehicle inertial parameters, anticipated disturbance environment, and 
the degree of response attenuation desired by gyric stiffning. An important end goal in designing 
an effective level of bias momentum is to mitigate the vehicle attitude response to unknown 
external disturbances. If the attitude response is measured in terms of mean-square response of 
the roll and pitch angular velocities in response to a specific bandlimited external disturbance 
torques, the approximate analytical formula given in equation 5.25 can be useful in terms of 
both sizing and qualitatively understanding the key factors and parameters. However, in any 
given physical application, it is the role of the design engineer to choose a physically suitable 
metric to quantify the vehicle attitude response and to secure a physically relevant model of 
the anticipated disturbance environment. 

Consider the following approximate but analytical MSR derived earlier in equation 5.25. 
Using the above expression, various special cases can be derived, including those given by 
equations 5.16, 5.21, 5.26, 5.27, and 5.28. In a typical sizing problem, one would expect to be 
given, at miniumum, the vehicle inertia, the disturbance bandwidth, and some idea of the level 
of damping in the system (c is either zero, “small” or “large”). With this basic information, the 
question is what level of bias momentum would be helpful in terms of mitigating the roll-pitch 
MSR to this assumed form of disturbances. 

To numerically examine the above sizing equations, consider the flying dual-spin system 
called the Dual Spin Test Device (DSTD) as shown in Figure 4. The DSTD is an indoor 



Figure 4: DSTD: A generic dual spin test device. 

flying platform at NASA Langley Research Center designed to validate dynamics and stability 
characteristics of a generic VTOL platform. It consists of a square honeycomb main body 
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with four ducted propellers fixed to the corners of the main body. A pair of control vanes 
are attached to the aft of each duct and a momentum wheel is attached to the bottom of the 
main platform. For more details on the DSTD, see [19], [20], [21], Table 1 shows the values for 
nominal parameters used for sizing a momentum wheel. 


Parameter 

Symbol 

Value 

vehicle mass 

771 

16.8 kg 

vehicle length, width, height 

h, hi h 

0.9, 0.9, 0.3 m 

vehicle roll, pitch, & yaw inertias 

hi hi h 

0.59, 0.58, 1.15 kg — m 2 

design disturbance torque variance 

e[r T r] 

14 N 2 m 2 

design disturbance torque bandwidth 

{BW) t 

3.2 Hz 

system damping 

C 

small but unknown 


Table 1: Nominal parameters values for DSTD. 


6.1 Dependence on System Damping 

Figure 5 shows the dependence of MSR on the levels of bias momentum and damping levels 
based on equation 5.25. Specifically, viscous damping coefficients of c = [0.01, 0.1, 1, 10, 100] 
are shown. The general trend is that the MSR decreases with increasing h B and damping 



Figure 5: MSR dependence on bias momentum for different levels of damping, 
in the system. Notice that based on the level of bias momentum, there appears to be three 
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different trends in the response. In Figure 5, the dashed-dotted vertical line indicates the 
angular momentum level such that 

h B = ( BW) t I -v=> x a — 1 

At this angular momentum level, the precession frequency A 0 is the same as the disturbance 
bandwidth (BW) T . This level can be viewed as a boundary between the “low” and the “inter- 
mediate/high” levels of bias momentum. 

For “low” bias momentum levels (on the left of the dash-dotted line), given a level of damping, 
the MSR is insensitive to different levels of bias momentum. In this case where x a > 1, the 
disturbance bandwidth is greater or equal to the precession frequency so that for small damping, 
the vehicle response will be dominated by a lightly damped precession mode. Once the lightly 
damped precession mode is excited and dominates the response, increasing the disturbance 
bandwidth or decreasing the bias momentum level will not change the MSR significantly. The 
special case responses given in equation 5.27 for small damping and equation 5.28 for large 
damping corresponds to this case. 

At intermediate or high momentum levels (on the right of the dash-dotted line) where x Q < 1, 
the disturbance bandwidth is less than the precession frequency so that as x a decreases, the 
precession mode will be excited even less, with a factor of 1/(1 — xl), as seen in equations 5.21 
and 5.26. Hence, at an intermediate h B level, where x Q ~ 1, the MSR is very sensitive (apparent 
even in log-log scale) to the level of bias momentum. Finally, at a high level of h B , the MSR 
decreases at two orders of magnitude per an order of magnitude increase in h B level. 


6.1.1 Validation of Approximate Analytical Formula 

To validate the approximations given in equation 5.25, these predictions are compared with 
detail simulation results from a fully nonlinear dynamical system response of the generic flying 
dual-spin system as described earlier. Based on a simulated sample disturbance torque of 
bandwidth 20 r/s and a damping level of c = 1, simulated responses at bias momentum levels 
of [1, 5, 10, 50, 500] Nms in the system are computed and represented by circles (o) in Figure 
5. Notice that the MSR based on the full non-linear simulation is close to the predicted MSR 
(represented by the dash-dot line), which validates the response approximation formula as given 
in equation 5.25. 

Figure 6 shows the angular velocity time responses in the roll axis for different levels of bias 
momentum in the vehicle system. Notice the different scales used to show the responses. It is 
clear that higher bias momentum levels can significantly attenuate angular velocity response of 
the vehicle. 

To simulate the errors incurred by violating the assumption of equal rotational inertias in 
the roll and pitch axes, the full nonlinear dynamics with J 2 //i = 0.7 is considered, and the 
results are shown in Figure 5 as square (□) symbols. Again, the simulation results, with the 
additional error due to unequal vehicle inertias, indicate a good match between the predicted 
response and the detailed simulation MSRs. 
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Figure 6: Roll angular velocity responses to disturbances at bias momentum levels of 

[5, 10, 50, 500] Nms and damping level of c = 1.0. 


Finally, notice from Figure 5 that the high degree of match in the MSR, especially at the 
higher bias momentum levels (or more specifically x Q -C 1), is due to the condition of bias 
momentum dominance as discussed earlier in equation 2.13. When the bias momentum level is 
gyroscopically dominant , the roll and pitch dynamical response becomes essentially linear, even 
at large angular velocities. 


6.2 Dependence on Vehicle Inertia 

Figure 7 shows the dependence of MSR on bias momentum levels for the following levels of 
vehicle inertia, [1, 2, 5, 10, 50]/ where / = 0.58 kg — m 2 . The damping coefficient of c = 1.0 
is assumed for all cases. As is clear from Figure 7, in general, the MSR value decreases with 
increasing bias momentum, and at higher bias momentum levels, the gyroscopic effects dominate 
and produces a second order decrease in MSR with increasing bias momentum. 

The dependence of MSR on the level of vehicle inertia, for a fixed level of bias momentum 
is slightly more complicated. At lower bias momentum levels, say h B = 5 Nms, the MSR 
decreases with increasing levels of vehicle inertia. However, at higher bias momentum, for 
example at h B = 20 Nms, it is evident from Figure 7 that the MSR for a vehicle with inertia 
/ is predicted to have a much smaller MSR than corresponding vehicles with larger inertias 
of 21, 5 1, and 10/. This counter-intuitive trend, that smaller inertia vehicles can give smaller 
responses than larger inertia vehicles for a fixed level of disturbances and bias momentum, 
follows from the fact that for a given level of bias momentum, smaller inertia vehicles more 
readily exhibit gyroscopic effects, as captured by the key parameter value of, x a < 1. In other 
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Figure 7: MSR dependence on bias momentum for different levels of vehicle inertia, at c = 1.0. 

words, the attenuation due to gyroscopic effects is greater than that due to reduction of motion 
due to an increase in inertia. For the extreme case with the very large inertia of 50/, the MSR 
finally becomes slightly lower than that for a vehicle with inertia /, at h B = 20 Nms. 

Notice that at h B = 20 Nms, only the smallest inertia vehicle with inertia / is expected 
to exhibit gyroscopic effects with a precession frequency of A 0 = h B /I = 34.5 r/s, which 
is significantly higher than the disturbance bandwidth, (BW) T = 20 r/s, and equivalently 
x 0 = 0.58 < 1. All the larger inertia vehicles, for example, a vehicle with inertia 10/ will have 
a much lower precession frequency of A 0 = h B /I = 3.4 r/s so that x Q = 5.9 > 1, will not exhibit 
significant attenuation due to gyroscopic effects. 

The dashed-dotted vertical line in Figure 7 indicates the angular momentum level such that 
(for a given vehicle inertia) the precession frequency A c is the same as the disturbance bandwidth 
(BW) t . This critical line increases with increasing vehicle inertia. This means that for larger 
inertia vehicles, higher levels of bias momentum are needed to induce sufficient gyric stability. 


6.3 Design Summary for DSTD 

Figure 8 shows a summary of relevant DSTD characteristics in the selection of the bias momen- 
tum level for disturbance rejection, based on Figure 5. A design range for the level of angular 
momentum needed is shown in Figure 8 along with the selected DSTD design of 17 Nms cor- 
responding to x D = 0.69 for a design disturbance as specified in Table 1. Its predicted response 
in terms of the standard deviation of the roll and pitch angular velocities is approximately 
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h_ (NMs) 

D 


Large damping effects dominate 
over gyroscopic and inertia effects 


Design Range 

Insensitive to damping levels 
Inverse square dependence; to BM 
Away from transient region 
Minimal BM level 


Gyroscopically dominant response 

£[ “ T ^mh® e[TTr] 


efurw] i 


2 cI(BW) Lightly damped precession response 


Figure 8: Summary of bias momentum sizing for DSTD. 

17.4 deg/ s. Of course the design response can be reduced further as necessary by implementing 
higher levels of angular momentum, for example, if the designed angular momentum is doubled 
to h B = 34 Nms, the predicted response drops to 6.7 deg/s. Notice also that the predicted 
response is directly proportional to the standard deviation of the torque disturbance levels and 
the design disturbance level of 3.7477m may be overly conservative. 

The salient features of the design range indicated in Figure 8: 

• The design range of h B lies in the region where the designed Precesssion frequency, A 0 is 
greater than the designed disturbance bandwidth, (BW) T , i.e. , x a < 1. 

• The lower end of the design range is slightly greater than h B — I ■ ( BW) r = 10.9 Nms, to 
avoid the large transient increase in response by the factor (1 — x%). The high end of the 
design range is driven by the need to limit the bias momentum level required, to mitigate 
weight and size requirements. 

• The design range lies in a region where the gyroscopic effects dominante the response and 
the response MSR has an inverse square dependence on the level of angular momentum. 

• In the design response range, the MSR is insensitive to damping levels in the system, with 
the exception that large damping effects dominate over gyroscopic and inertia torques, a 
physically unlikely scenario. This is an important advantage because damping levels for 
systems are typically difficult to model and predict accurately. 
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A physical realization of a bias momentum wheel designed for 17 Nms is shown in Figure 9 and 



Figure 9: A momentum wheel fabricated for DSTD. 

consists of a flywheel having a diameter of 0.39 m, ring mass of 1.08 kg (6.4 % of gross weight 
of DSTD vehicle), and spinning nominally at 4000 rpm. As much as technologically feasible, 
the mass of the flywheel can be reduced by either increasing the diameter of the flywheel or 
the flywheel spin rate. 

In sizing bias momentum for the purpose of mitigating a vehicle response to external dis- 
turbances, the simplest approach is to either increase damping sufficiently (Figure 5), and/or 
increase inertia sufficiently (Figure 7) such as to satisfy a certain desired level of MSR for the 
particular vehicle. In both approaches, dual-spin stabilization using bias momentum is not 
necessary. However, increasing vehicle damping or inertia sufficiently for a free-flying vehicle 
may be detrimental to performance or even be physically unrealizable, for a VTOL application. 

The alternative to the above predicament is to “roll-off” the MSR with sufficiently large bias 
momentum, i.e., dual-spin stabilize the vehicle. A generally desirable goal then is to select a 
bias momentum level which is on the 1/h 2 slope, i.e., for a given vehicle inertia and disturbance 
bandwidth, choose a bias momentum level such that the precession mode frequency will exceed 
the disturbance bandwidth to some chosen margin, or based on equation 3.8, 

h B > (1 + margin) \Zl 1 I 2 (BW) T (6.1) 

Of course, avoiding exciting the precession frequency will be particularly important if the pre- 
cession motion is undamped or only very lightly damped. By going further to the right side of 
the vertical line (where x a = 1), the MSR becomes less and less dependent on inertia and damp- 
ing changes. In the extreme case where x Q -C 1, the above simplifying behavior approaches 
static approximations, meaning that the disturbance bandwidth is well below the precession 
frequency. In general, one could size the bias momentum based on a requirement to stay below 
a target mean square response as given in equation 5.25 or below a target frequency response 
as given in equation 4.14. 
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7 Simulation Results 


This subsection show various results based on DSTD vehicle simulation that demonstrates the 
advantages in using bias momentum in a flying platform. The significant improvements in 
the vehicle attitude control robustness to disturbances are not limited to open loop hovering 
conditions. The simulation results indicate that improvements in performance can extend to 
closed loop response, command tracking, and even payload changes during flight. Interestingly, 
simulation results also demonstrates that the use of bias momentum is most dramatic under 
larger disturbances in payload or winds. 

This subsection describes simulation results based on full nonlinear models that are con- 
structed from first principles, bench test data refined model parameters, sampled and quantized 
measurements, and saturation limits on actuators. In particular, the dynamical models for the 
ducted propcller/vane, whose development are based on bench tests, includes asymmetrical ef- 
fects which are difficult to visualize due to their 3-dimensional characteristics. This MATLAB- 
Simulink simulation model also incorporates models of servo dynamics that regulate fan speeds 
and vane angles. 

In the simulation cases that follow, the vane angles and the fan speeds are actively used for 
stabilization and control which reflects the current state-of-the art approach based on vector and 
differential thrusting. Initially, the simulation focused on the effects of bias momentum in the 
stabilization and control performance of the flying test platform. Gyroscopic stability effects 
are highlighted by comparing a system with significant level of bias momentum (BM cases) 
with a corresponding system with no bias momentum (NBM). A difference in this study is the 
independent control of all four propeller speeds and their respective vane angles, as opposed to 
constraining it to work in pairs and collectively, as in conventional helicopter control [22], [23] 
or practically all VTOL flying platforms [24]. 


7.1 Crosswind and turbulence model 

To simulate wind disturbances during hovering or low airspeed operations, crosswinds are as- 
sumed to be impinging on a vertical cylinder attached to the upper part of the platform, as 
shown in Figure 10. The resulting aerodynamic forces on this cylindrical column is used to 
simulate the crosswind and turbulence on the vehicle. The empirical model given in [25] for 
aerodynamic forces over a cylindrical column is used. The resulting disturbance model consists 
of forces and torques in the pitch and roll axes. Two sets of winds are simulated, namely, 
“light” winds, and “strong” winds. The “light” winds simulate “calm to gentle breeze” with 
sample mean wind speeds of 7 Knots and standard deviation of 5. However for “strong” winds, 
intended to simulate “calm to strong breeze/near gale”, the simulated sample mean wind speed 
is 15 Knots with a standard deviation of 18, as seen in Figure 11. 
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Figure 10: Simulated turbulent crosswinds on DSTD. 
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Figure 11: Simulated strong wind. 
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7.2 Controller Law Synthesis 

Figure 12 shows a schematic of the model used for control law design and simulation. The 


Vane Vane angle 



trim command 


Figure 12: DSTD model used for control law design and simulation. 

DSTD mathematical model with fan and servo models has 24 states consisting of 3 linear 
momenta, 3 angular momenta, 3 attitudes, 3 positions, 8 fan motor states, and 4 vane servo 
motor states. The 8 fan motor states include for each ducted fan, first order PI controllers to 
regulate individual fan speeds which are assumed governed by first order dynamics. Based on 
LQG control theory (see for example [16]), feedback controllers are designed to stabilize and hold 
the platform attitude and also track the translational velocity and angular velocity commands. 
In all cases, 4 vane angles and 4 fan speeds were commanded. For control during both hovering 
and velocity tracking, 12 states including 3 attitudes, 3 angular velocities, 3 positions, and 
3 linear velocities were used for feedback by the LQG controller. In this simulated design 
study, feedback controllers for the BM and NBM cases are synthesized based on the linearized 
model about hovering trim, which is identical to the trim conditions at a constant translational 
velocity if aerodynamics due to vehicle airspeed are not included, as assumed in this study due 
to their low dynamic pressures. For a fair comparison between BM and NBM cases, the same 
LQG weights were used in the controller design 

A key difference in the BM and NBM cases is that in the BM case, each pair of fans 
are counter rotated to mitigate both inertial and aero-resistance fan drag torques due to fan 
rotations, as compared to same-direction rotation of all four fans in the BM case. Consequently, 
the trim conditions in BM and NBM cases are necessarily different since the net aero-resistance 
fan drag torques are different, and must be statically canceled by vane-deflected forces. 

Figure 13 shows the tradeoff between performance cost and control effort for a parameterized 
set of LQG controllers undergoing gust response during hover over a period of time Tf ina i. The 
quadratic cost function used in the LQG control law included six states for position and velocity 
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Figure 13: Control effort (abscissa) vs Performance cost (ordinate) for a parameterized set of 
LQG controllers undergoing gust response during hover. 


and six states for attitude and angular rates. The different controllers were parameterized by 
scaling the quadratic control penalty term consisting of four vane angles and four fan speeds 
excursions about their individual trim values. For ease of comparison of closed loop performance 
among different LQG control laws and for both cases with BM and NBM, consider only a subset 
of variables. Specifically, the performance cost in terms of attitude ($) and position (£j) errors 
are 


m = / £[(a - A) 2 + a - iifidt 


(7.1) 


i — 1 


and the control effort is defined by a cumulative control activity (CCA) function as follows 


CCA(t) 



(7.2) 


The above measure of control effort is a time integral of the sum of squares of all four control vane 
angular rates (6Q. The abscissa and ordinate in Figure 13 represents CCA(Tf ina i ) and J(Tf ina {) 
respectively. In addition to quantifying the level of control activity, the above control effort is 
intended to capture the level of control vane activity, which is seen as the primary source of 
control surface induced unsteady aerodynamics due to its thrust vectoring role. Albeit limited, 
this tradeoff study indicates that controllers are significantly more effective for a system endowed 
with significant levels of bias momentum (asterisk line) in terms of actuator effort necessary 
for a given level of performance cost. Although it is clear from this trade study that the NBM 
case requires significantly higher levels of vane control activity, the effect of controller induced 
unsteady aerodynamics is not included due to a lack of reliable model. However, one would 
expect to see a larger degradation of performance for the NBM cases. 


28 


7.3 Hovering 

7.3.1 Open Loop Response to Light Winds. 

The open-loop attitude responses of the initially trimmed system for the two cases are shown in 
Figure 14. In the NBM case (red line), the vehicle rapidly loses orientation and goes unstable, 
but in the BM case (blue line), the vehicle holds attitude inspite of a small gyric roll and pitch 
angle response to the light winds. However, the platform drifts off slowly without feedback 
compensation due to a steady wind component in the light winds. 
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Figure 14: Open loop response to light winds, 
Case (red line). 
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7.3.2 Closed Loop Response to Light Winds. 

The closed-loop responses to light winds during hover for both cases are shown in Figure 
15. Although both cases show stable response, the angular response for the BM case (blue 
line) is clearly less oscillatory and therefore superior to the NBM case (red line). Apparently, 
under these simulated light wind conditions, the angular response amplitudes for the NBM case 
appears to be small enough to be tolerated during an actual flight. The time history of the 
cumulative control effort (bottom right of Figure 15), shows that the case with BM actually 
requires less vane actuation than in the corresponding NBM case. As expected, the vane and 
fan actuator response did not saturate for both cases, under light winds. 
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7.3.3 Closed Loop Response to Strong Winds. 


The closed-loop responses to strong winds during hover for both cases are shown in Figure 16. 
The controller for NBM case was judiciously scaled to mitigate actuator saturation. Both cases 
remains stable but the angular response for the BM case (blue line) is clearly less oscillatory 
and therefore superior to the NBM case (red line). However, the significantly larger angular 
response amplitudes for the NBM case will likely be unacceptable or even dangerous during an 
actual flight. The time history of the cumulative control effort shows that the case with BM 
requires less vane actuation and it gives far superior attitude hold than the corresponding NBM 
case. 

To examine what happens if the controller for NBM case is not judiciously scaled to miti- 
gate actuator saturation, resulting possibly to a more agressive feedback control, an additional 
simulation was done with the new control law. It turns out that the close loop system loses 
control after a few seconds. More specifically figures 17 and 18 shows the time histories of the 
vane angle and fan speed responses respectively. The vane and fan actuators saturate resulting 
in loss of vehicle under strong winds. 





Figure 15: Closed loop response to light winds. BM Case (blue line), NBM Case (red line). 
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Figure 16: Closed loop response to strong winds. BM Case (blue line), NBM Case (red line). 




Figure 17: Closed loop response of control vanes to strong winds, more agressive control. BM 
Case (blue line), NBM Case (red line). 
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Figure 18: Closed loop response of fan speeds to strong winds, more agressive control. BM 

Case (blue line), NBM Case (red line). 
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7.4 Command Tracking 


This section outlines a sample of simulation results of command tracking under strong winds, 
specifically translational velocity commands while holding a trim attitude. It should be noted 
that the demonstrated inherent advantage of bias momentum in endowing directional stability 
to the platform does not necessarily improve its maneuverability. Hence, the simulation results 
in this subsection is intended to examine this issue quantitatively. 

For the following simulations, LQG based controllers are designed to track a trajectory 
path defined by a commanded translational velocities about trim with zero angular velocity. 
The input commands are shaped to mitigate the discontinuity in a step command. Figure 19 
shows the response to translational velocity tracking commands (dashed green line) for both 
cases subject to strong winds. The accompanying excursions in the angular velocities (and 
hence the vehicle attitude) are shown in Figure 20. It is seen that for both cases, the vehicle 
can successfully track translational velocity commands but BM case (blue line) gives superior 
performance. 



Figure 19: Vehicle translational velocity response during trajectory tracking. BM Case (blue 
line), NBM Case (red line), Command (dashed green line. 


The corresponding vane and fan actuator time histories during tracking are shown in Figures 
21 and 22, respectively. The case with BM (blue line) clearly needs less vane and fan control 
activity but as observed earlier gives smaller velocity tracking errors. 
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Figure 20: Vehicle angular velocity response during trajectory tracking. BM Case (blue line), 
NBM Case (red line). 





Figure 21: Vane angle response during tracking. BM Case (blue line), NBM Case (red line). 
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Figure 22: 
(red line). 
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Fan speed response during command tracking. BM Case (blue line), NBM Case 
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7.5 Robustness to payload variations 

7.5.1 Simulated Mass and Inertia changes. 

To simulate variable payloads, assume a step increase of 20 % (3.36 kg) of the gross weight of the 
platform. A spherical uniformly distributed mass is assumed dropped at 0.15 m away from the 
center of mass location on the platform at 5 sec and is assumed to be dropped off (or removed) 
from the platform at 20 sec. The inertia matrix and total mass of the system is significantly 
changed as a result of the off-centered loading. This off-CG centered loading simulation is of 
particular interest due to its difficulty and its practical significance for a hovering platform 
where passengers or payloads can be loaded and unloaded with maximum flexibility without 
compromising vehicle stability while hovering. In this part of simulation, disturbances and 
measurement errors are not included to focus entirely on the effects of these payload changes. 


7.5.2 Closed loop response during hover. 

Figures 23 and 24 show the position and attitude responses to mass and inertia changes for 
both cases. It is seen that the platform position and attitude excursions are significantly better 
in the BM case as shown by the blue line. 
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Figure 23: Position response to step changes in payload. BM Case (blue line), NBM Case (red 
line) . 
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Figure 24: Attitude response to step changes in payload. BM Case (bine line), NBM Case 

(red line). 


Significant differences are also noted in their corresponding vane and fan actuator histories as 
shown in Figures 25 and 26 respectively. It is seen that the fan speeds and vane angles oscillate 
significantly and rapidly in the NBM case resulting in the platform oscillations shown previously. 
Unfortunately, the actual consequences on flight performance due to this apparently more 
unsteady control vane activity for the NBM case, remains to be verified during an actual flight 
because a reliable mathematical model is not available for this class of unsteady aerodynamics. 
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Figure 25: Control vane responses to step changes in payload. BM Case (blue line), NBM 

Case (red line). 
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Figure 26: Fan speed responses to step changes in payload. BM Case (blue line), NBM Case 
(red line). 
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7.6 Summary of Simulation Results 


Simulation results demonstrate the potential performance of a new stability and control tech- 
nology, based on trying to mitigate control surface actuator induced unsteady aerodynamics 
effects in feedback control. Its ultimate goal is to enable statically unstable, disturbance sen- 
sitive, thrust levitated, low-speed or hovering vehicles to operate more safely than is possible 
today, in turbulent conditions and under large payload variations. Simulation results indicate 
that generic VTOL platforms can be endowed with significant directional stability with a proper 
amount of bias momentum, similar to a dual-spin stabilized spacecraft but under a significantly 
different operating environment. The results also indicate that this open loop directional sta- 
bility carries over to closed loop if the controllers are designed properly. Specifically, it was 
found that vehicle attitude control robustness can be significantly improved (over a no-bias- 
momentum vehicle) during hovering, particularly under strong turbulent winds or significant 
payload variations. In addition, command tracking performance at low airspeeds can also be 
improved significantly particularly under turbulent winds, despite the presence of somewhat 
complicated gyroscopic coupling. These performance improvements are due to control laws 
that are based on proper modeling of the system dynamics. These predicted improvements in 
robustness obtained from detailed simulation studies are expected to be more evident during 
actual flight test, with the physical presence of control induced unsteady aerodynamics. 


8 Conclusions 


A new methodical approach to sizing bias momentum for dual spin hovering vehicles has been 
proposed which is based on the roll and pitch responses to disturbances. The approach is 
heuristic in nature and addresses a very practical but fundamental issue of bias momentum 
sizing for passive gyroscopic stability and the influence of the major factors including the as- 
sumed disturbance environment, nominal vehicle rotational inertia, and the designed level of 
attitude response of the vehicle to disturbances. Further work on realizing hardware system 
for storing a large amount of bias angular momentum is recommended. Clearly more detailed 
tradeoff and optimization in the design of momentum wheels tailored for air vehicles are also 
needed to improve performance further. Although spin stabilization is an established technol- 
ogy in spacecraft attitude control, flight validation is strongly recommended due to the vastly 
different operational environment faced by air vehicles. 
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Appendix 


Indefinite integrals used: 
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Taylor series of tan 1 x is 
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